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Abstract 

The results of a detailed histogram Monte-Carlo study of critical-fluctuation 
effects on the magnetic-field temperature phase diagram associated with the 
hexagonal Heisenberg antiferromagnet with weak axial anisotropy are re- 
ported. The multiphase point where three lines of continuous transitions 
merge at the spin-flop boundary exhibits a structure consistent with scaling 
theory but without the usual umbilicus as found in the case of a bicritical 
point. 
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The frustration of antiferromagnetically coupled sites in a triangular array gives rise to 
noncolinear spin order and magnetic-field temperature phase diagrams which exhibit a rich 
variety of structure. 0i In the case of a hexagonal lattice with weak (c-axis) axial anisotropy, 
a novel type of multicritical point occurs at which three lines of continuous transitions merge 
at the spin-flop and paramagnetic phase boundaries. Although the full phase diagram has 
been observed only in the quasi-one-dimensional hexagonal antiferromagnet CsiVzC^i'S 
similar behavior is expected in0 RbNiCh, CsMnI 3 , and CsiVz.Br3.lll Gross features of the 
experimental results for CsNiCl^ have been reproduced by the analysis of a phenomeno- 
logical Landau-type free energy constructed from symmetry arguments.! In zero field, two 
continuous transitions occur, identified by the components of the spin vector M.§ As the 
temperature is lowered from the paramagnetic phase 1, a linear phase 2 with M z (where 
z||c) is stabilized at TVi, and at Tjv2 < TJvi an elliptically polarized phase 3 with M z and 
M x occurs. In a magnetic field applied along the c axis, a first-order spin-flop transition 
to a helically polarized phase 4 with M x and M y is found. Each of the ordered states is 
characterized by a period-three modulation in the basal plane and a period-two structure 
along the c axis. (The helical phase 4 is similar to the well-known 120° spin struture of 
triangular antiferromagnets). Scaling analysis of the multicritical point where phase 1, 2, 
3, and 4 meet suggests that the behavior of all the three lines of continuous transitions 
are governed by the same crossover exponent 0.0 It is the purpose of the present work to 
examine in detail this behavior by means of accurate histogram Monte Carlo simulations of 
the anisotropic Heisenberg antiferromagnet. 

The easy-axis Hamiltonian studied here is given by0 

H = J,, £ S, ■ Sj (1) 

<ij> 

where J\\ > and J± > are the nearest-neighbor antiferromagnetic exchange interactions 
along the c axis and in the basal plane, respectively, D < is the single-ion anisotropy, and 
H is magnetic field applied along the c axis. Although Jn >> J± for CsNiCl^, we consider 



here the isotropic case J» = J± = 1 for simplicity and since Monte-Carlo simulations of 
models with strongly anisotropic parameters require considerably more computing effort. 
For this reason, we also chose a value for D = —0.2, which is small enough for our pur- 
poses; the ground states which occur as a function of H are consistent with those observed 
experiment ally.0 (In contrast with the case D = — 1, for example, which yields a phase dia- 
gram with a completely different structure. 00) The present work serves to compliment and 
extend the Monte Carlo studies performed at H=0 on this and similar Hamiltonians,00 as 
well as one cursory examination of the phase diagram.0 

This study was made in an attempt to accurately estimate the phase-boundary lines 
close to the multicritical point at (T m , H m ). It was anticipated to be a numerically challeng- 
ing problem since fluctuations involving all three components of M are important in this 
region of the phase diagram. The Ferrenberg-Swendsen histogram method of Monte Carlo 
simulations offers the possibility of the precise determination of transition points by the 
temperature at which extrema occur in thermodynamic functions for finite-size systems.0 
Relevant components of the staggered susceptibility, defined according to the components of 
the order parameter M involved in the transition of interest ,0 were used for this purpose. 
Simulations were performed on a lattice of size 12 x 12 x 12. Runs of 1.2 x 10 6 Monte Carlo 
steps per spin were made, with the initial 2 x 10 5 steps discarded for thermalization. For a 
given value of magnetic field, single histograms were made at one or more T to ensure that 
the maxima in the susceptibility occured close to at least one simulation temperature. 

The results shown in Fig. 1 confirm the general stucture determined by the phenomeno- 
logical Landau model! as well as by a molecular-field treatment of the Hamiltonian (1).0 (In 
the latter case, two different types of linear and elliptical states (2A, 2B, 3A, and 3B) were 
distinguished by a relative phase angle, as in the Monte Carlo studies at H=O.00 Such 
distinctions are beyond, and not relevant to, the goals of the present work.) The detailed 
behavior near the multicritical point at T m = 0.915(5), H m = 2.62(4), however, is quite 
different from both mean-field results. There is clear indication that both the 1-2 and 1-4 
transition lines approach this point with slopes that are the same in magnitude as that of 



the 3-4 spin-flop line. This is an effect of critical fluctuations and is predicted by scaling 
theory.0 Less clear is that the 2-3 transition line exhibits the same predicted tendency. A 
remarkable feature of the 1-4 line is its initial curvature to the left as the field increases. 
In the usual case of the bicritical point associated with unfrustrated antiferromagnets with 
weak axial anisotropy, only two critical lines are involved.0 These also approach this point 
asymptotically with slopes that are the same in magnitude as the spin-flop line and are 
governed by the same crossover exponent (p. In contrast with the present case, both of these 
lines have always been found to approach the bicritical point from the right, forming an 
umbilicus structure.0 This behavior has also been observed in Monte Carlo simulations!!! 
There appears to be no argument from scaling theory, however, that relates the signs of 
the initial slopes of critical lines emanating from a multicritical point. We note that this 
opposite-slope behavior has not been observed in CsNiCl^. This is not surprising since 
the model parameters used in here are not relevant for quasi-one-dimensional materials, a 
feature which may obscure this unusual effect. 

In conclusion, these Monte Carlo results demonstrate that significant critical-fluctuation 
effects are associated with this novel multicritcial point. While there are strong symmetry 
arguments to support the conclusion made in Ref. [TO] that the 1-2 and 2-3 transitions belong 
to the xy and Ising universality classes, respectively, the nature of the 1-4 transition, and 
that of the multicritical point itself, remain somewhat unsetteled. The interpretation made 



in Ref. |10| is that these transitions are related to the new chiral universality classes proposed 
by Kawamura.0 (The resulting crossover-exponent value ~ 1.04 is not inconsistent with 
the results of Fig. 1.) Azaria et a&~§! have argued that such transitions exhibit nonuniversal 
critical behavior, where a first-order, mean-field tricritical or 0(4) universality can occur. It 
is not clear what type of scaling behavior for the multicritcal point can be expected in these 
cases. 
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FIGURES 

FIG. 1. Phase diagram with H||z||c near the multicritical point at T m = 0.915(5),iT m = 2.62(4) 
as determined by Monte Carlo simulations (points). Indicated are the paramagnetic phase 1, linear 
phase 2 with M z , elliptical phase 3 with M z , M x and spin-flopped helical phase 4 with M x , M y . 
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